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ABSTRACT 


—  ^  The  purpose  of  this  thesis  was  to  determine  whether  or 
not  there  are  phenomenological  differences  between  dynamic 
derivatives  calculated  from  F-15A  rotary  balance  data  and 
data  from  other  sources. 

To  make  this  determination ,  two  additional  data  sets 
were  obtained:  (1)  the  F-15A  design  phase  stability  deri¬ 
vative  data,  and  (2)  the  F-15A  production  phase  stability 
derivative  data.  The  lateral  dynamic  derivatives  were  then 
compared  through  derivatives  of  the  lateral  moments  with 
respect  to  the  rotation  rate  about  the  velocity  vector  (wind 
vector) . 

The  conclusion  of  the  project  was  that  differences  ex¬ 
ist  between  the  data  sets,  but  that  the  dominant  character¬ 
istics  were  the  same  for  all  of  the  data  sets;  the  differ¬ 
ences  in  the  data  were  not  indicative  of  basic  (phenomeno¬ 
logical)  differences  in  the  data  itself.  Therefore,  the 
contention  that  oscillatory  rates  affect  determination  of 
the  dynamic  derivatives  was  not  substantiated  by  this  study. 


THE  EFFECT  OF  CONSTANT  VS  OSCILLATORY  RATES 


ON  DYNAMIC  STABILITY  DERIVATIVES 


I.  INTRODUCTION 


Historical  Survey 

Traditionally,  aircraft  stability  and  contro1  is 
determined  from  a  linearized  system  of  differential 
equations  of  motion.  The  linearization  is  obtained  by 
expanding  the  force  and  moment  terms  (as  functions  of  the 
motion  and  control  variables)  in  truncated  Taylor  series 
about  an  equilibrium  point.  The  partial  derivatives  of  the 
force  and  mcment  coefficients  with  respect  to  the  motion 
variables  are  called  aerodynamic  stability  derivatives  and 
have  been  commonly  used  over  the  years  to  predict  aircraft 
stability  characteristics. 

The  stability  derivatives  are  often  further  subdivided 
into  "static"  and  "dynamic"  stability  derivatives.  The 
static  derivatives  are  those  taken  with  respect  to  V,  <*  , 
and  s  .  All  other  derivatives  are  considered  "dynamic". 
Basically,  the  static  derivatives  are  relatively  "easy  and 
inexpensive"  to  determine  from  wind  tunnel  tests  while  dy¬ 
namic  derivatives  are  relatively  "difficult  and  expensive" 


to  determine 


Because  of  the  difficulty  and  expense  involved  in  the 
extensive  dynamic  testing  required  to  obtain  the  dynamic 
derivatives  empirically  (a  requirement  necessary  for  their 
accurate  determination  over  a  broad  range  of  flight  condi¬ 
tions)  ,  the  aviation  industry  has  shied  away  from  doing 
this  testing.  The  industry  uses  flight  test  data  to  deter- 
mine  them  instead.  Prior  to  the  mid-1930's,  a  combination 
of  the  inaccuracy  of  then-current  test  equipment,  the 
relatively  low  cost  of  full  scale  aircraft  prototypes,  and 
the  relatively  conventional  designs  produced  encouraged  the 
designers  to  estimate,  or  in  some  cases  ignore,  the  dynamic 
derivatives  during  the  design  phase.  Stability  problems 
attributable  to  the  inaccuracies  in  their  estimates  were 
corrected  through  modifications  during  or  after  the  pro¬ 
duction  phase. 

With  the  advent  of  the  more  sophisticated,  high  per¬ 
formance  designs  introduced  as  a  result  of  the  Second  World 
War,  increased  emphasis  was  placed  on  predicting  the  sta¬ 
bility  and  performance  characteristics  of  new  aircraft 
prior  to  production  (or  even  first  flight) .  But  the  cost 
of  dynamic  testing,  was  still  prohibitive.  In  response  to 
the  need  for  data,  a  myraid  of  new  estimation  techniques 
was  developed  and  put  into  use  by  the  industry  (1). 

The  techniques  developed  for  VWII  prevailed  through 
the  following  decade.  In  the  mid-1950's,  two  new  factors 


emerged  to  discourage  the  use  of  dynamic  testing  to  deter¬ 
mine  the  dynamic  derivatives  in  the  design  phase:  (1) 
analysis  of  air  operations  in  the  Korean  War  indicated  tre¬ 
mendous  U.S.  supremacy  over  U.S.S.R.  aircraft  as  demonstra¬ 
ted  by  the  12-1  kill  ratio  achieved  in  air  to  air  combat 
(2:  110) ,  and  (2)  the  increasingly  rapid  development  of 

electro-mechanical  automatic  control  devices.  The  combina¬ 
tion  of  these  two  factors  made  it  easier  and  cheaper,  and 
therefore  more  attractive,  to  modify  new  designs  during  (or 
after)  production,  if  necessary,  to  make  them  meet  required 
stability  and  performance  criteria. 

In  the  last  decade,  a  combination  of  several  fac¬ 
tors  has  reawakened  interest  in  dynamic  testing:  (1) 
analysis  of  air  operations  in  Viet  Nam  showed  a  narrowing 
of  U.S.  supremacy  as  indicated  by  the  2.5-1  kill  ratio 
achieved  (2:  110) ,  (2)  the  reduction  in  excess  capability 

available  in  the  automatic  contol  systems  to  correct  inher¬ 
ent  design  problems  (due  in  large  part  to  the  expanded  use 
of  flexible  structures  throughout  modern  aircraft) ,  and  (3) 
the  development  of  high  speed/low  cost  computers  which  pro¬ 
vide  a  relatively  low  cost,  easy  way  to  acquire  and  reduce 
the  large  amounts  of  data  required  by  dynamic  testing. 

Without  the  technological  superiority  enjoyed  in  the 
past,  mistakes  in  the  design  phase  can  no  longer  be  cor¬ 
rected  in  the  production  phase  and  still  achieve  acceptable 


performance.  Therefore,  new  means  of  acquiring  the  infor 
mation  necessary  in  the  design  phase  in  a  cost  effective 
manner  must  be  evaluated. 


Determining  the  Dynamic  Derivatives 

There  is  a  variety  of  techniques  currently  in  use  to 
determine  the  dynamic  derivatives.  These  techniques  range 
from  the  relatively  simple  and  inexpensive  estimation  meth¬ 
ods  contained  in  the  DATCOM  (3)  and  similar  sources  (1,  4) 
to  the  very  complex  and  expensive  experimental  procedures 
developed  by  various  research  facilities  throughout  the 
world  (5) . 

Over  the  last  decade,  the  NASA-Langley  Research  Cen¬ 
ter  has  been  developing  a  technique  for  predicting  spins 
using  a  rotary  balance  installed  in  the  NASA-Langley  Ver¬ 
tical  Wind  Tunnel.  Birhle  (6),  in  outlining  the  history  of 
spin  related  research,  points  out  the  difficulty  past  in¬ 
vestigators  have  had  in  predicting  spins  and  the  trajec¬ 
tories  of  spinning  aircraft.  He  goes  on  to  say  that  the 
reason  the  former  investigators  were  unsuccessful  was  that 
the  data  necessary  to  compute  the  forces  and  moments  due  to 
the  steady  state  components  of  the  body  angular  rates  -  p, 
q,  and  r  -  must  come  from  a  constant  rotation  dynamic  test 
device,  like  the  rotary  balance,  which  was  not  formerly 
available.  For  accurate  trajectory  calculation,  he  con¬ 
tends,  the  rotational  velocities  must  be  broken  up  into 


steady  state  and  oscillatory  components  and  the  contribu¬ 
tions  of  each  to  the  forces  and  moments  determined  by  tests 
in  steady  state  facilities  (e.g.,  rotary  balance)  and  fa¬ 
cilities  that  are  not  solely  steady  state,  respectively. 


Purpose  of  This  Thesis 

If  it  is  true  that  steady  state  and  oscillatory  com¬ 
ponents  of  the  body  angular  rates  contribute  to  the  forces 
and  moments,  other  than  in  proportion  to  their  magnitudes 
and  directions,  then  the  dynamic  derivatives  computed  from 
steady  state  sources  and  sources  other  than  solely  steady 
state  will  be  different.  Further,  these  differences  must 
be  able  to  be  explained  due  to  phenomenological  differences 
in  the  flows  of  the  source  data  testing  apparatus.  The 
purpose  of  this  thesis  is  to  examine  dynamic  derivatives 
computed  from  a  data  set  acquired  from  a  solely  steady 
state  source  and  some  data  sets  acquired  from  non  solely 
steady  state  sources  to  determine  any  significant  differ¬ 
ences,  and,  if  they  exist,  the  phencmological  reasons  for 
these  differences. 


II.  Topic  Discussion 


This  section  introduces  and  discusses  topics  necessary 
for  the  analysis  carried  out  in  subsequent  sections.  The 
first  subsection  looks  at  the  characteristics  of  aerodynamic 
test  data  to  determine  guidelines  for  comparing  data  sets. 
Next,  the  phenomena  that  affect  the  forces  and  moments  due 
to  oscillatory  motion  are  examined  and  their  probable  con¬ 
tributions  due  to  oscillatory  rotation  rates  analyzed.  The 
nature  of  rotary  balance  force  and  moment  coefficient  data 
is  examined  in  the  third  subsection.  And,  finally,  the  dy¬ 
namic  derivatives  taken  with  respect  to  the  body  angular 
rates  are  analyzed. 

Aerodynamic  Test  Data 

Collection  and  analysis  of  aerodynamic  test  data  is 
necessary  to  experimentally  determine  the  values  of  aero¬ 
dynamic  parameters.  A  persistent  problem  with  this  pro¬ 
cess,  however,  has  been  that  the  data  collected  often  shows 
significant  scatter.  This  scatter  shows  in  a  limited  manner 
for  successive  trials  of  a  single  test  and  more  dramatically 
for  tests  conducted  under  a  variety  of  conditions  or  for 
dissimilar  tests  designed  to  determine  a  certain  parameter. 
As  a  result  of  the  scatter,  then,  the  existence  of  differ¬ 
ences  between  parameters  calculated  from  different  data 
sets  has  come  to  be  expected. 


When  the  objective  changes  from  determining  the  values 
of  aerodynamic  parameters  to  canparing  data  sets  through 
parameters  calculated  from  these  data  sets,  some  different 
implications  are  drawn  fran  the  scatter.  Simple  differ¬ 
ences  in  the  parameters  may  represent  scatter  rather  than 
differences  in  the  underlying  data  and  therefore  are  not 
sufficient  to  demonstrate  significant  differences  in  the 
data  sets.  To  demonstrate  such  differences,  variations 
between  the  major  characteristics  of  the  parameters  must  be 
shown. 

There  are  many  characteristics  of  such  parameters  that 
can  be  usefully  compared.  Among  them  are  slopes  and  magni¬ 
tudes,  points  of  change  of  slope,  and  consistency  (or  lin¬ 
earity).  The  choice  of  the  most  appropriate  characteristic 
to  use  is  dependent  on  the  particular  parameter  under  con¬ 
sideration. 

Once  the  characteristics  to  be  compared  have  been 
chosen,  demonstration  of  similarity  or  dissimilarity  is 
based  on  two  criteria:  (1)  the  parameters  compared  must 
predominantly  exhibit  the  major  characteristics,  and  (2) 
the  similarities  or  dissimilarities  identified  must  be  con¬ 
sistent  throughout  the  data  sets.  The  key,  then,  to  estab¬ 
lishing  similarity  or  dissimilarity  rests  with  showing  con¬ 
sistency.  Without  consistency,  demonstrations  of  similar¬ 
ity  or  dissimilarity  become  mere  demonstrations  of  scatter. 


Oscillatory  Rate  Phenomena 

There  are  two  phenomena  which  may  affect  the  aero¬ 
dynamic  forces  and  moments  as  a  result  of  oscillatory 
rotation  rates.  These  phenomena  are  unsteady  aerodyna¬ 
mics  and  aeroelasticity. 

Unsteady  aerodynamic  effects  result  from  the  delays 
in  the  response  of  the  flow  field  about  an  aircraft  to 
motions  of  the  aircraft.  When  the  aircraft  motions  are 
rapid  compared  to  the  response  of  the  flow  field,  tran¬ 
sient  forces  and  moments  are  applied  to  the  aircraft  while 
the  flow  field  adjusts  to  the  changing  conditions.  These 
forces  and  moments  caul  be  very  pronounced  when  the  motion 
is  due  to  the  structural  modes  of  the  aircraft. 

These  same  structural  modes  can  cause  large  forces  and 
moments  due  to  changes  in  the  incidence  of  the  wing.  When 
triggered  by  unsteady  forces  of  the  proper  frequency  on  the 
aircraft  surfaces,  the  surfaces  vibrate,  sometimes  with  re¬ 
latively  large  amplitude  deflections.  These  deflections  can 
cause  changes  in  the  incidence  of  the  wing  leading  to  addi¬ 
tional  forces.  (If  these  additional  forces  are  of  the  pro¬ 
per  frequency,  resonance  results  in  very  large  amplitude 
deflections.) 

The  body  angular  rates  are  rigid  body  modes  of  the  air¬ 


craft.  The  characteristic  frequencies  are  on  the  order  of 
6-8  rad/sec  (7).  The  characteristic  frequencies  for  un- 


steady  aerodynamic  phenomena  and  aeroelastic  phenomena  are 
in  the  range  20  rad/sec  and  up  (8).  Since  both  of  these 
phenomena  require  exitation  frequencies  of  the  same  order 
as  their  own  frequencies,  they  will  not  be  excited  by  the 
body  angular  rates.  As  a  result,  these  phenomena  will  not 
contribute  to  the  forces  and  moments. 

Rotary  Balance  Force  and  Moment  Data 

The  rotary  balance  used  at  the  NASA-Langley  Research 
Caiter  provides  for  an  aircraft  model  to  be  mounted  in  the 
spin  tunnel  at  any  desired  angles  of  attack  and  sideslip 
and  then  spun  about  an  axis  parallel  to  the  velocity  vector 
while  the  forces  and  moments  are  measured  and  recorded  (6) . 
In  normal  mountings,  the  bank  angle  is  zero.  The  angle  of 
pitch  is  preset  so  it  is  unaffected  by  the  spin;  thus  the 
pitching  velocity'  is  zero.  This  means  that  the  spinning 
velocity,  and  therefore  the  force  and  moment  data,  is  a 
function  of  the  rolling  and  yawing  body  rates,  p  and  r, 
only. 

The  spinning  velocity  of  the  aircraft  is  the  overall 
rotation  rate,  n.  In  general,  n  is  a  combination  of  the 
roll  and  yaw  rates,  p  and  r,  which  cannot  be  independently 
controlled.  Therefore,  the  resulting  force  and  manent  data 
can  only  be  given  as  a  function  of  this  combination.  It 
is  not  possible  to  use  the  spin  tunnel  force  and  moment 
data  to  find  the  rates  of  change  (partial  derivatives)  of 


the  forces  and  moments  with  respect  to  the  roll  and  yaw 
rates  individually. 

The  Dynamic  Derivatives 

As  stated  before,  the  dynamic  derivatives  are  all  of 
the  stability  derivatives  except  those  taken  with  respect 
to  V,  «,  and  fi.  However,  only  those  derivatives  taken  with 
respect  to  the  body  angular  rates  -  p,  q,  and  r  -  are  ap¬ 
plicable  to  this  report. 

The  combination  of  the  six  force  and  moment  coeffi¬ 
cients  and  the  three  angular  body  rates  has  the  potential 
to  produce  18  derivatives.  The  uncoupling  of  the  longitu¬ 
dinal  and  lateral  modes  reduces  this  by  half  to  nine,  how¬ 
ever,  and  the  fact  that  the  rotary  balance  data  is  not  a 
function  of  the  pitching  rates  reduces  this  number  still 
further  to  six. 

The  yawing  moment  coefficient  is  normally  larger  in 
magnitude  than  and  proportional  to  the  side  force  coeffi¬ 
cient.  Therefore,  the  information  provided  by  the  two  co¬ 
efficients.  is  redundant.  Since  the  lateral  modes  are  nor¬ 
mally  more  dependent  on  the  moments  than  the  forces,  this 
report  will  not  examine  the  side  force  derivatives.  This 
leaves  the  four  moment  coefficient  derivatives. 

Finally,  since  these  four  derivatives  cannot  be  inde¬ 
pendently  derived  from  rotary  balance  data,  only  the  com¬ 
bined  derivatives  (the  derivatives  of  the  moment  coeffi- 
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cients  with  respect  to  the  overall  rotation  rate)  can  be 
used  in  this  analysis. 
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III.  APPROACH 


To  accomplish  the  thesis'  objective,  three  things 
needed  to  be  done:  (1)  choose  data  sets,  (2)  choose  the 
comparisons  to  make  between  the  data  sets,  and  (3)  make  the 
comparisons.  This  section  deals  with  the  first  two  activ¬ 
ities. 

Data  Set  Selection 

The  F-15  was  chosen  as  the  aircraft  to  use  as  the  test 
case  because  it  is  the  only  operational  Air  Force  aircraft 
to  have  been  tested  on  the  rotary  balance  (10) .  An  opera¬ 
tional  aircraft  was  desired  because  a  larger  data  base 
exists  for  operational  aircraft  than  for  most  other 
aircraft. 

Two  data  sets  were  chosen  for  ccmparison  with  the  ro¬ 
tary  balance  data.  The  first  set  (11)  was  one  for  which  the 
dynamic  derivatives  were  estimated  (12  -  14) .  The  second 
set  (15)  was  derived  from  production  phase  flight  test  data. 
(Note:  It  would  have  been  desirable  to  have  a  data  set  de¬ 
rived  from  oscillatory  wind  tunnel  testing  -  to  be  a  source 
of  solely  oscillatory  data  -  but  because  the  Air  Force  pol¬ 
icy  for  the  last  three  decades  has  been  to  forego  this  type 
of  testing,  the  data  is  not  available.) 

Choice  of  Conparisons 

The  rotary  balance  data  set  provides  force  and  moment 


coefficient  data  as  a  function  of  rotation  rate  (from  which 
stability  derivative  parameters  can  be  computed)  for  one 
Mach/Reynolds  number  combination  and  a  wide  range  of  angles 
of  attack.  The  design  phase  and  production  phase  data  sets 
provide  the  stability  derivatives  for  a  large  number  of 
Mach/Reynolds  number  combinations  but  for  a  limited  angle 
of  attack  range. 

Reynolds  number  (Re)  is  determined  by  the  relation: 

Re  *  flVl/*i  (1) 

which,  when  written  in  terms  of  Mach  number,  becomes: 

Re  =  ,*Mal/M  (2) 

The  quantity  a  %  decreases  with  increasing  altitude,  so 
for  a  given  Mach  number,  Reynolds  number  decreases  with  in¬ 
creases  in  altitude  (see  Table  I). 

The  stability  derivative  data  in  the  design  and  pro¬ 
duction  phase  data  sets  are  either  independent  of  altitude 
(Reynolds  number)  or  tend  toward  a  limiting  value  with  in¬ 
creasing  altitude  (decreasing  Reynolds  number) .  For  Mach 
numbers  of  0.6  and  below  this  limiting  value  is  very  nearly 
reached  by  80,000  ft.  Therefore,  this  altitude  value  was 
used  to  compute  Reynolds  numbers  for  the  design  phase  rigid 
data  following  the  example  of  the  production  phase  data 
(but  note  was  taken  that  the  Reynolds  number  for  the  rigid 
data  can  be  an  arbitrarily  lower  value;  in  the  limit  it  is 
zero) .  A  summary  of  Reynolds  numbers  for  the  data  sets  is 
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given  in  Table  II. 

The  most  meaningful  comparisons  are  those  for  compara¬ 
ble  Mach  numbers,  Reynolds  numbers,  and  angles  of  attack. 
For  this  analysis,  the  comparisons  were  limited  to  data  for 
incompressible  Mach  numbers  (0.6  and  less),  Reynolds  num¬ 
bers  within  approximately  one  order  of  magnitude,  and  the 
common  angle  of  attack  range. 

TABLE  II 

Data  Set  Mach  and  Reynolds  Numbers 


Report 

Mach  Number 

Reynolds  Number 

NASA  CR  3478 

0.022 

0.211  X  106 

MDC  A0502 

0.2 

0.908  X  106 

0.6 

2.725  X  106 

MDC  A4172 

0.3 

1.363  X  106 

0.6 

2.725  X  106 

IV.  DATA  ANALYSIS 


The  Data  Sets 

The  three  data  sets  were  produced  over  a  period  of  ap¬ 
proximately  twelve  years  (roughly  1969-1981).  As  a  result, 
the  reports  represent  tests  of  slightly  different  configura¬ 
tions  of  the  F-15A.  (Some  of  the  more  noticeable  changes 
were  the  addition  of  the  leading  edge  snag  on  the  horizon¬ 
tal  stabilizer  and  downsizing  of  the  speed  brake  which  oc¬ 
curred  in  the  early  seventies.)  In  addition,  some  of  the 
tests  configured  the  aircraft  differently,  for  example,  in- 
including  or  not  including  air  to  air  missiles.  The  effects 
of  these  variations  in  the  configuration  are  small,  however, 
having  a  negligible  affect  on  this  analysis. 

Rotary  Balance  Data  (Appendix  A) .  The  rotary  balance 
data  is  raw  moment  coefficient  data  plotted  as  a  function 
of  angle  of  attack,  non-dimensional  spin  rate,  and  spin 
radius.  The  portion  of  the  data  package  used  for  the  com¬ 
parison  was  limited  to  the  data  for  the  basic  configura¬ 
tion  with  no  controls  deflected  and  for  the  angle  of  at¬ 
tack  range  8°-  35°  which  has  a  spin  radius  equivalent  to  six 
feet  on  the  full  scale  aircraft.  (Note:  The  data  for  an¬ 
gles  of  attack  greater  than  35°  was  taken  with  the  spin 
radius  set  to  zero.) 

Design  Phase  Data  (Appendix  B) .  The  design  phase  data 
consists  of  plots  of  the  stability  derivatives  as  a  function 


of  angle  of  attack,  Mach  number,  and  altitude.  (Note: 

Only  the  data  for  a  rigid  aircraft  was  used.) 

Low  angle  of  attack  (less  than  approximately  26 8 )  and 
high  angle  of  attack  (approximately  26°  to  35°)  data  are 
presented  separately  and  there  are  some  inconsistencies  in 
the  data  in  the  interface  range.  These  inconsistencies, 
while  noticeable,  do  not  appreciably  affect  the  results. 

The  data  used  for  comparison  was  limited  to  the  two 
Mach  numbers  0.2  and  0.6  and  the  angle  of  attack  range  8” 
to  35°. 

Production  Phase  Data  (Appendix  C) .  The  production 
phase  data  is  presented  in  basically  the  same  format  as  the 
design  phase  data.  The  data  is  not  presented  differently 
for  two  angle  of  attack  ranges,  so  it  doesn't  suffer  fran 
the  inconsistencies  of  the  other  data  set. 

The  data  used  for  comparison  was  limited  to  the  two 
subsonic  Mach  numbers  0.3  and  0.6,  the  angle  of  attack 
range  8°  to  35°,  and  80,000  ft  altitude. 

Computing  the  Comparison  Parameters 

As  stated  earlier,  the  rotary  balance  force  and  moment 
coefficient  data  is  a  function  of  a  combination  of  the  roll 
rate,  p,  and  the  yaw  rate,  r.  As  a  result,  the  individual 
stability  derivatives  with  respect  to  these  rates  cannot, 
in  general,  be  determined  from  the  rotary  balance  data, 


while  the  derivatives  with  respect  to  the  overall  rotation 


rate  can.  The  appropriate  parameters  for  comparison,  then 
are  the  derivatives  with  respect  to: 


Rotary  Balance  Comparison  Parameters.  The  comparison 
parameters  were  calculated  from  the  rotary  balance  moment 
coefficient  data  in  Appendix  A  for  the  entire  angle  of  at¬ 
tack  range  (8 8  -  90°)  using  a  linear  least  squares  curve 
fit. 

Design  and  Production  Phase  Comparison  Parameters. 

These  stability  derivatives  were  calculated  from  the  sta¬ 
bility  derivative  data  in  Appendices  B  and  C  using  the 
relation: 

Cj  =  C_-  cos  «  +  sin«  i  =  l,n  (3) 

1n  *r 

which  is  developed  in  Appendix  D. 

Presentation  of  the  Comparison  Parameters 

The  computations  of  the  comparison  parameters  along 
with  the  results  were  tabulated  and  the  results  plotted  (see 
Appendix  E).  Plots  for  the  angle  of  attack  range  8°  to  35° 
were  made  for  each  data  set/Mach  number /parameter  combina¬ 
tion;  in  addition,  a  plot  of  each  parameter  was  made  for  the 
rotary  balance  data  for  the  entire  angle  of  attack  range. 

The  results  for  each  parameter  for  all  of  the  data  set/Mach 
Number  combinations  for  the  angle  of  attack  range  8°  to  35° 
are  summarized  in  Figures  1  and  2. 


V.  Results 


This  section  analyzes  the  plots  discussed  in  Section 

IV.  It  looks  at  the  plots  of  C.  analyzing  the  similar i- 

in 

ties  and  differences  of  the  rotary  balance,  design  phase, 

and  production  phase  data.  It  then  analyzes  the  plots  of 

Cn  in  the  same  manner. 
nn 

Rolling  Moment  Derivative  (C.  )  Plots.  The  overall  trends 

_ 

of  the  plots  of  C1  in  Figure  1  show  that  the  magnitude  of 

n 

of  the  derivative  increases  from  a  large  negative  value  at 
8°  angle  of  attack  to  a  slightly  negative  value  at  25°  an¬ 
gle  of  attack  and  then  stays  relatively  constant  from  25° 
to  35°.  This  dominant  trend  is  seen  for  all  the  plots. 

The  rotary  balance  data  taken  with  a  six  foot  spin 
radius  does  not  show  the  bend  at  25°,  but  shows  it  delayed 
to  30°.  However,  the  data  taken  with  zero  spin  radius 
(Figure  49)  shows  the  bend  at  the  expected  25°. 

The  design  phase  data  shows  the  trend  change  at  25° 
as  noticeably  as  the  rotary  balance  data  but  does  not  show 
the  remarkable  linearity  from  8°  to  25°  that  the  rotary  bal¬ 
ance  data  does.  Since  C,  is  predominantly  C,  for  low  an- 

1n  XP 

gles  of  attack,  the  nature  of  C-i  should  account  for  the 

difference.  The  plots  of  C.  (Figures  17  and  23)  show  that 

1P 

it  does  not  change  uniformly  with  angle  of  attack  but  shows 
a  sharp  increase  from  10°  to  15°  for  the  0.2M  data  and  from 


20°  to  25°  for  the  0.6M  data.  This  step  function  nature  is 

not  reflected  in  the  rotary  balance  data;  the  step  shifts  to 

lower  angle  of  attack  ranges  for  lower  Mach/Reynolds  number 

combinations  and  may  have  shifted  out  of  the  comparison 

range  for  the  data  in  this  report. 

The  production  phase  data  once  again  change  in  trend 

at  25°  angle  of  attack.  The  range  from  8 °  to  25 °,  however, 

is  characterized  by  peaks  and  valleys  in  stark  contrast  to 

both  the  rotary  balance  data  and  the  design  phase  data. 

Both  C  and  C.  (Figures  31  -  34)  are  highly  nonlinear  and 
iP  xr 

the  nonlinearity  of  C.  again  dominates  C.  for  both  0.3M 

xp  xn 

and  0.6m  although  the  contribution  of  C.  is  readily  detect- 

Ir 

able  from  8°  to  12°  for  0.3M. 

Yawing  Moment  Derivative  (C_  )  Plots.  Following  the  pattern 

_  nQ 

seen  with  C.  ,  the  plots  of  Cn  in  Figure  2  show  a  change 

in  trend  at  approximately  25°.  While  the  trend  following  25° 

is  not  as  consistent  or  obvious  for  Cn  as  it  was  for  C,  , 

n 

it  is  a  pronounced  change,  inmediately  seen  in  all  but  one 
data  set. 

Paralleling  the  data  for  C.  ,  the  rotary  balance  Cn 

xn  « 

data  for  the- range  8°  to  25°  is  surprisingly  linear.  It 
follows  a  negative  trend  from  8 0  to  20°,  bottoms  out  from 
20°  -  25°,  and  increases  slightly  from  25°  -  35°. 

The  design  phase  data  show  a  peak  between  8°  and  25° 
at  15°  for  0.2M  and  at  12°  for  0.6M  before  picking  up  the 


O-  Rotary  Balance  Data 
A-  Design  Phase  Data  0.2M 
V-  Design  Phase  Data  0.6M 
□  -  Production  Phase  Data  0.3M 
O"  Production  Phase  Data  0.611 


the  downward  trend  of  the  rotary  balance  data.  Both  bottom 
out  between  20 3  and  25°  before  the  0.2M  data  picks  up  a  neu¬ 
tral  trend  and  the  0.6M  data  picks  up  a  negative  trend.  The 
breed?  at  25°  is  not  as  pronounced  for  the  0.6M  data  but  is 
noticeable.  The  trends  after  25°  for  both  0.2M  and  0.6M 
differ  frcm  that  for  the  rotary  balance  data. 

The  peak  between  8°  and  25°  for  both  sets  of  data  is  a 
result  of  following  Cn^,  whose  behavior  from  0°  to  16° 
(Figures  19  and  25)  is  not  reflected  in  the  rotary  balance 
data.  The  production  phase  data  show  the  same  trends  as  the 
design  phase  data  for  the  range  8  °  to  25  °  but  the  downward 
trend  from  12°  to  20°  is  much  more  pronounced  (it  plunges). 
As  a  result  of  being  so  much  more  negative  at  25°  than  the 
other  data  sets,  the  trend  from  25°  to  35°  is  strongly  posi¬ 
tive  in  order  to  return  to  the  range  near  zero  of  the  other 
data  sets. 

Overall,  the  outstanding  characteristic  of  all  the 
data  is  that  the  data  shows  one  trend  from  8°  to  25°  and 
another  trend  after  25°.  This  occurs  regardless  of  the 
Mach  and  Reynolds  number  of  the  data  set  and  the  parameter 
looked  at.  The  dominant  trend  prior  to  25’  is  positive  for 
and  negative  for  Cn  ;  after  25®,  and  Cn  both  tend 
toward  a  nearly  constant  value. 


VI .  Conclusions 


The  analysis  has  shown  that  there  are  both  similar¬ 
ities  and  differences  between  the  three  data  sets.  The  dom¬ 
inant  trends  are  exhibited  by  all  three  data  sets  but  the 
characteristics  of  these  trends  differ  within  each  set  to 
a  greater  or  lesser  degree. 

The  analysis  does  not  show  a  phenomenological  differ¬ 
ence  between  the  data  taken  on  the  rotary  balance  and  the 
other  data  sets.  Differences  exist,  but  the  differences 
are  not  consistent;  the  rotary  balance  data  for  C  for  8°  - 
25°  is  very  similar  to  the  design  phase  data  but  very 
different  for  the  production  data  for  this  range  while  all 

three  are  similar  for  the  range  from  25°  to  35°.  For  Cn  it 

n 

is  just  the  opposite. 

In  surtmary,  significant  differences  in  the  lateral  dy¬ 
namic  derivatives  were  neither  predicted  by  the  theory  nor 
seen  in  the  actual  test  results.  Therefore,  this  report 
concludes  that  the  contention  that  there  is  a  significant 
difference  between  data  taken  for  steady  state  rotation 
rates  and  oscillatory  rotation  rates  is  not  substantiated 
by  this  analysis. 
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Rotary  Balance  Data:  Cn  vs  nb/2V  for  SR 


Rotary  Balance  Data:  Cn  vs  nb/2V  for  SR 
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Figure  17.  Design  Phase  Data:  C,  vs*  for  0.2M 


Figure  1ft.  Design  Fhase  Data:  C.  vs  «*•  for 


Figure  19.  Design  Phase  Data:  C  vs  *  for  0.2M 
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Phase  Data: 


Figure  22.  Design  Phase  Data:  C_,  vs  «•  for  0.611 


Figure  23.  Design  Phase  Data:  C..  vs  for  0.6M 


Figure  24.  Design  Phase  Data:  C1  vs  for  0.6M 


Figure  26.  Design  Phase  Data:  C  vs  *  for  0.6ff 


Figure  27.  Production  Phase  Data:  C  vs  *-  for  0.3M 
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Figure  29.  Production  Phase  Data:  C  vs  «■  for  0.3M 
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Figure  33.  Production  Phase  Delta: 
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Figure  35.  Production  Phase  Data:  C.  vs  a  for  0.311 
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Figure  37.  Production  Phase  Data:  C  vs  a  for  0.3M 


Production  fhase  Data:  C  vs  a  for  0.6M 
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APPENDIX  D:  Derivation  of  the  Relation  for  Cj 

a 

When  a  model  is  tested  on  the  rotary  balance,  it  is 
mounted  initially  with  the  x  axis  of  the  body  axis  system 
aligned  with  the  velocity  vector  (i.e.,  opposite  the  wind 
vector) .  Since  this  corresponds  to  an  aircraft  pointing 
along  the  gravity  vector,  the  orientation  of  the  y  and  z 
axes  is  arbitrary.  The  angle  of  attack  (a)  is  defined  as 
the  angle  between  the  x  axis  and  the  wind  vector,  and  so  it 
is  zero  in  this  initial  mounting. 

For  testing,  the  model  is  set  at  an  angle  of  attack 
and  then  rotated  about  an  axis  parallel  to  the  wind  vector. 
This  rotation,  when  viewed  in  the  body  frame,  is  composed 
of  a  roll  rate  and  a  yaw  rate.  To  find  an  expression  for 
these  rates  in  terms  of  the  overall  rotation  rate, n ,  it 
is  necessary  to  do  a  1-2  (9,0)  axes  rotation  from  the  wind 
(inertial)  axes  frame  to  the  body  axis  frame  (  0  here  is  the 
rotation  about  the  x,  axis  so:  n). 

Designating  the  w^  to  be  the  wind  axes  frame  unit  vec- 

A 

tors,  ij  to  be  the  intermediate  axes  frame  unit  vectors, 

A 

and  bj  to  be  the  body  axes  frame  unit  vectors,  the  trans¬ 
formations  become: 


D-l 


ss 


cos  a  0  -sin 


sin  a  ■  0  cos  a 


The  overall  rotation  in  the  wind  axes  frame  is: 


_*b/i 

<3  =  aw, 


Using  the  transformation  relation: 


PI  -  w-M 


the  rotation  in  the  intermediate  frame  is: 


-+b/i  _  ? 

oj  3  Al. 


Transforming  to  the  body  frame: 


m  • 


o  o  r 

n  L  J,- 


so  the  rotation  in  the  body  frame  is: 


K  A  ,  A 

=  cos «  bt  +  sm«  b3 


Since  the  body  rates  are  defined  by: 


p  b,  +  q  b2  +  r  b3 


the  body  rates  beccme: 


WW5WK5WB 


r  =  nsin» 

or  in  non-dimensional  form: 


(12) 


pb/2V  »  (ab/2V)  cos  « 


(13) 


rb/2V  -  (Ab/2V)sin« 

The  effect  on  the  force  and  manent  coefficients  of  the 
overall  rotation  rate  must  be  the  same  as  the  sum  of  the 
effects  of  the  body  rates  so: 


q  />'  =  C.  p’  +  C.  r'  (14) 

*  xp  1z 

where  a' ,  p' ,  and  r'  are  nondimens ional  and: 


C,  =  aCi/  a(kb/2 V)  (15) 

xk 

i  *  Nf  Af  1^  m,  n 
k  38  p' ,  r' ,  a' 

Substituting  for  p'  and  r': 

Cj  £  -  C.  a' cos  *+  C.  A'sin*  (16) 

XP  h 


and  dividing  through  by  a': 


Appendix  E:  Tabulations  and  Plots  for  Each 
Parameter/Mach  Number/t)ata  Set 
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Rotary  Balance  Data: 
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Table  VII 


Design  Phase  Data:  Computation  of  C,  vs  *  (8°  -  35°) 
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Table  VIII 


Design  Phase  Data:  Computation  of  Cn  vs  *(8°  -  35°) 
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Table  XI 


Production  Phase  Data:  Computation  of  C,  vs  *  (8°  -  35°) 
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Production  Phase  Data:  Computation  of  Cn  vs  *  (8®  -  35°) 
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